Cs CPT magnetometer for cardio-signal detection in unshielded environment 
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We present first, encouraging results obtained with an experimental apparatus based on Coherent 
Population Trapping and aimed at detecting biological (cardiac) magnetic field in magnetically 
compensated, but unshielded volume. The work includes magnetic-field and magnetic-field-gradient 
compensation and uses differential detection for cancellation of (common mode) magnetic noise. 
Synchronous data acquisition with a reference (electro-cardiographic or pulse-oximetric) signal allows 
for improving the S/N in an off-line averaging. The set-up has the relevant advantages of working 
at room temperature with a small-size head, and of allowing for fast adjustments of the dc bias 
magnetic field, which results in making the sensor suitable for detecting the bio-magnetic signal at 
any orientation with respect to the heart axis and in any position around the patient chest, which 
is not the case with other kinds of magnetometers. 



I. INTRODUCTION 

Spectroscopy of transitions involving long-lived levels 
in alkali atoms has shown impressive potentialities in the 
field of high resolution and high sensitivity magnetom- 
etry since the birth of coherent spectroscopy [3, 0, 
In the last years, many research groups, with different 
experimental techniques like non linear magneto-optical 
rotation spin-exchange relaxation free - Faraday rota- 
tion 0], double resonance optical pumping and Co- 
herent Population Trapping (CPT) 0,11] have obtained 
important records in optical magnetometry sensitivity, so 
to prove that atomic magnetometers are presently com- 
petitive even with SQUID magnetometers in terms of 
sensitivity, and can find applications in fields like geo- 
magnetism monitoring, testing of materials, testing of 
fundamental symmetries of Physics and bio-magnetism 
detection. 

The first mapping of the cardiomagnctic field with 
a Rubidium double resonance magnetometer has been 
demonstrated in ref. Q, and good agreement with the 
magnetocardiogram obtained by a SQUID magnetocar- 
diograph has been found [Io|. 

The unique experiment of M. Romalis and coworkers 
(Tl| . recently showed that laser magnetometers can per- 
form records in registration of such weak (few hundreds 
fT) magnetic fields as those produced by the human brain 
activity. 

A fundamental request for actually reaching the sen- 
sitivity level needed for bio-signal detection is to reduce 
dramatically the environmental magnetic field. Usually 
high magnetic permeability materials like /i-metal are 
employed to build small volumes inside of which it is 
possible to have very effective shielding of the external, 
undesired magnetic fields. 

The high cost of such magnetically clean chambers rep- 
resents a strong limitation for the diffusion of high sensi- 
tive magnetometry in clinical application, also because of 



their delicateness and possible magnetization with con- 
sequent requirement of periodical demagnetizing treat- 
ments. 

The main limits of working in unshielded environmen- 
tal conditions, typical in a scientific laboratory, are given 
by the presence of strong magnetic field gradients, ac 
magnetic fields in the volume of sensor cell and by the 
magnetic noise produced by other human activities and 
ionospheric phenomena. 

Field gradients affect the sensitivity limit of the instru- 
ment by introducing a broadening of the detected reso- 
nance line. Time-dependent magnetic fields represent a 
further additional source limiting the magnetic field de- 
tection sensitivity. 

The possibility of achieving high sensitivity, simply 
by passively quenching magnetic noise (placing thick Al 
plates around the magnetic sensor), has been shown in 
ref.[T2|. There a sensitivity of the order of lpT/\/Hz is 
obtained thanks to a sophisticated feedback system that, 
keeping around zero the total magnetic field value in the 
volume of a K cell (heated to 170°C), suppresses the 
broadening due to spin-exchange-collisions. 

In the present work, we show that a cheap system for 
eliminating spatial inhomogeneities of the background 
magnetic field, joined to an accurately balanced differen- 
tial detection scheme, allows for detecting the magnetic 
signal produced by the heart activity without the use of 
neither expensive and bulky shielding chambers nor alu- 
minum shields that attenuate high-frequency magnetic 
noise. The device consists in a full optical sensor operat- 
ing at room temperature, placed inside a frame of com- 
pensation coils that can be easily accessed by a human 
body. 
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II. BASICS OF THE MEASUREMENT 
PRINCIPLE 

Magnetic field optical detection is performed by creat- 
ing CPT[1,[i^] on Zeeman sublevels of the F=3 hyperfine 
ground state of Cs atoms. 

This is an all optical technique as it does not require 
the presence of coils in the proximity of the sensor cell 
in order to produce direct RF magnetic excitation. This 
feature makes easy to optimize the sensitivity to any ori- 
entation of the magnetic field, thus allowing for measur- 
ing any component of it. CPT is produced when circu- 
larly polarized laser radiation, resonant with a D2 optical 
transition, is frequency modulated exactly at the Zeeman 
frequency splitting of the ground Zeeman sublevels. The 
CPT signal can be seen as a resonant decrement of the 
absorption in the optical transition. The measure of the 
resonant frequency u>l gives directly an absolute mea- 
surement of the modulus of the magnetic field interacting 
with the atomic probe. 

This kind of magnetic field measurement is absolute 
and self-calibrated because the magnetic field strength 
B is given by: B — , where is the Bohr magne- 
ton, gp is the Lande factor of the considered ground-state 
and uil is the resonance (Larmor) frequency. The magne- 
tometer is a scalar one and the estimation of the magnetic 
field modulus is obtained by scaling the Larmor frequency 
by atomic constants, and only minor deviations may oc- 
cur [Til ], introducing small systematic errors. Long time 
living coherence between Zeeman sublevels allows for cre- 
ating very narrow (few Hz) resonances and then for get- 
ting very sensitive magnetic field estimations. Cs optical 
transitions can be excited by means of diode laser radia- 
tion, that can be easily modulated in frequency by means 
of the modulation of its injection current. Magnetic field 
intensities in the geophysical range, correspond to fre- 
quency modulations of the order of 100 kHz. This gives 
the possibility of performing high precision magnetic field 
measurement while operating in the presence of the lo- 
cal Earth magnetic field, by employing RF modulation 
of the laser current [si fTBL ITq] . 

The presented experimental investigation is dedicated 
to the detection of a small, time varying field in the pres- 
ence of about 10 6 times stronger dc bias field, for any 
relative orientation of the two fields. As any scalar mag- 
netometer, our setup is sensitive to the component of the 
varying magnetic field in the direction of the bias field 
(see also Section HVA|) . 



III. SENSOR SET-UP 



tical isolator, the temperature stability system (Peltier 
junction and the thermistor), and the fiber collimator is 
contained into a compact butterfly housing of 40 cm 3 . 
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FIG. 1: Experimental setup. PMOF: polarization main- 
taining optical fiber, CL:collimating lens, IBS: intensity beam 
splitter, NF: neutral filters, PD: photodiode. 

The beam exiting from the fiber is collimated then split 
by an intensity beam splitter (IBS1). The beam reflected 
by IBS1 passes trough a vacuum reference cell and the 
transmitted light is detected in order to monitor the laser 
frequency tuning with respect to the Doppler line. The 
beam transmitted trough IBS1 is split again by a sec- 
ond intensity beam splitter IBS2. The beam reflected 
by IBS2 is made parallel to the transmitted one after 
reflection on a mirror M and plays the role of a second 
arm of the sensor. The two arms interaxial separation 
A is chosen to be about 7 cm. The polarization of both 
the beams is adjusted by means of two A/4 plates. Both 
beams cross, in different points, a single cell, containing 
Cs and 5 Torr of At as a buffer gas. The use of a sin- 
gle cylindrical cell, 9 cm in length and 1.5 cm in radius 
guarantees equal conditions for CPT creation (equal Cs 
and Ar vapor density in the two laser-atoms interaction 
volumes). All the component represented in FigfT] are 
mounted on a non magnetic plate and the overall dimen- 
sions of the sensor are 25 x 12 x 4 cm 3 . 

In order to reduce power broadening, laser power is 
reduced to less than 50 fiW in each sensor arm using a 
set of neutral filters. A non magnetic translating blade is 
placed after the cell, in front of one (the most sensitive) 
of the two photo-diodes, in order to have the possibility 
of fine balancing the detected photo-currents. 

In order to increase the signal to noise ratio of the 
detected signal, phase sensitive detection is used. At this 
aim, we impose a 20 kHz frequency modulation on the RF 
frequency modulating the laser current and then detect 
the component of the signal oscillating at this frequency. 



The sensor set-up is sketched in FigfTJ Laser source 
is a single-mode edge-emitting pigtail laser (A= 852 nm) 
with 15 mW of laser power and an intrinsic line-width 
of less than 5 MHz. Laser light is coupled to a single- 
mode polarization-maintaining fiber 10 m in length and 
the laser head, containing the laser chip, a 40 dB op- 



IV. SENSITIVITY OPTIMIZATION FOR 
CARDIO-SIGNAL DETECTION 

The magnetic field produced by a human heart is char- 
acterized by a maximum intensity of about 100 pT in the 
close vicinity of the chest, and the main features are well 



3 



reconstructed if sampled in a bandwidth of at least 30 Hz 

Magnetic heart-beat spatially confined distribution, of 
the order of the heart dimensions, allows for performing 
differential detection in a rather small volume in front of 
the patient's chest, while the possibility of triggering the 
acquisition on other signal, taken for example from the 
electro-cardiogram (or even from a simple pulse oxime- 
ter), allows for effective noise rejection in an off-line anal- 
ysis. 

As introduced before, unshielded environmental con- 
ditions limit the instrumental sensitivity essentially be- 
cause of the temporal and spatial fluctuation of the back- 
ground magnetic field. Field gradients in our laboratory 
have the typical magnitude of about 100 nT/cm depend- 
ing on the presence of magnetic material both in the in- 
strumentation and in the building structure. Alternating 
magnetic fields contribute instead, mainly with 50 Hz and 
higher odd order harmonics, with a typical rms amplitude 
of the order of 100 nT depending on the vicinity of elec- 
tric lines and power supply transformers. For frequencies 
higher than the lock-in cut-off frequency, ac fields deter- 
mine an effective broadening of the CPT line, while for 
lower frequencies they directly introduce time fluctuation 
of the CPT resonance center. Optimization of the mag- 
netometer sensitivity is fundamentally performed on the 
basis of the spatial and temporal characteristics of the 
bio-signal to be measured. 

In order to improve the common mode noise rejection it 
is important to have, in absence of the magnetic source 
to be measured, exactly the same resonance profile at 
the two inputs of the differential lock-in amplifier. As 
shown in Ref.[7j small differences in laser power, polar- 
ization, beam-size contribute significantly to the unbal- 
ance by changing the amplitude and the shape (mainly) 
of the CPT resonance. Difference between the magnetic 
field values inside the two laser-atom interaction volumes 
shifts instead the positions of the resonances centers. 

The cardiac signal measurement procedure consists in, 
first, adjustment of amplitude and width of the CPT sig- 
nals in the two arms separately (single input mode of the 
lock-in amplifier) and, second, fine optimization of those 
parameters by looking at the maximum common mode 
noise cancellation in the differential signal (differential 
input mode of the lock- in amplifier). 



A. Field inhomogeneities compensation 

As mentioned in Section [TTJ measuring the modulus of 
B in the presence of a strong dc component, leads to de- 
tect only the variations AB of the field in the component 



different components of the measured field and eventu- 
ally reconstructing the entire vectorial signal. 

The typical configuration for a cardio-signal detection 
measurement is sketched in Fig[2] In this particular con- 
figuration both x and y components of the laboratory 
magnetic field are compensated (the residual magnetic 
field is in the 0.1 f/T range) and thus, the z component 
of the human heart magnetic beat is being detected. 



parallel to the dc, bias field since A 



B ^rfy-AB. 

By zeroing two of the three spatial component of the 
local background magnetic field, one is able to measure 
one selected component of the signal produced by the bi- 
ological source. This gives the possibility of studying the 



FIG. 2: Schematics of the apparatus. The selected component 
of the background bias field Bo, parallel to z in the scheme, 
determines the direction of the detected component of the bio- 
signal, while laser beams separated by A are directed along 
x. Dipoles are used for gradients compensation. 

This detection scheme has important advantages in 
view of simplifying the compensation of magnetic field 
inhomogeneities. As explained in deeper detail in the 
Appendix, in fact, only 3 currents for the compensation 
of the spatial gradients inside the cell volume, are needed 
in addiction to the 2 currents employed for the zeroing of 
two of the three component of the laboratory magnetic 
field. These three currents flow respectively in one anti- 
Helmholtz pair and in two small coils playing the role of 
magnetic dipoles. 

After the optimization of the gradient compensation 
we get a strong reduction of the line broadening due to 
magnetic field inhomogeneities. Fig. [3] shows the CPT 
resonance signal registered with and without the com- 
pensation of spatial gradients. Without compensation 
a best fit analysis gives a linewidth (FWHM) of about 
1.7 kHz which is reduced down to less than 300 Hz in 
case of gradients compensation. 



B. Differential measurement procedure and noise 
characterization 



Depending on the spatial configuration of the local 
magnetic field in correspondence of the two light-atoms 
interaction volumes, we are able to cancel different noise 
contributions at the cost of increasing the uncorrelated 
electronic noise in the detection stage by namely a fac- 
tor of V2- Other noise sources, like laser frequency noise 
and amplitude noise are systematically canceled from the 
signal provided that the two arms are balanced. If, fur- 
thermore, both arms of the sensor give CPT resonance 
exactly at the same resonant frequency, then even the 
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FIG. 3: CPT line-shapes produced with circular polarization 
and magnetic field orientation orthogonal to the laser beam 
(along z direction). Circles: static gradients are not com- 
pensated. Fitted line-width is 1.67 kHz. Squares: optimized 
compensation of the gradient components. Fitted line- width 
is reduced to 280 Hz. 



common mode magnetic field noise, generated for exam- 
ple by all far magnetic sources is eliminated. 

In Fig. [4] is presented the magnetometric noise spec- 
tral density for different operating conditions of the de- 
vice. Trace (a) is relative to a single arm operation: laser 
noise and magnetic common mode noise contribute en- 
tirely. Trace (b) is obtained in differential input acqui- 
sition, where the two CPT resonances occur at frequen- 
cies separated by more than the CPT linewidth. In this 
case one mainly gets the cancellation of the laser fre- 
quency and intensity noise. Evident gain in the reduc- 
tion of residual noise is obtained in the differential acqui- 
sition when CPT resonances detected by the two sensors 
are accurately overlapped (trace (c)). Trace (d) and (e) 
are recorded respectively in differential and single input 
mode, when the laser frequency is tuned out of the sin- 
gle photon (Doppler) resonance. It is evident that the 
uncorrelated electronic noise increases in the differential 
detection. 



V. CARDIO SIGNAL DATA ANALYSIS 

The signal-to-noise (S/N) ratio is not high enough to 
make single cardiac pulse directly observable, so that off- 
line analysis must include both linear filtration and av- 
eraging of the data. We use a standard technique for 
acquiring a reference signal simultaneously with the mag- 
netometric data. Our digital lock-in amplifier (Stanford 
SR830) allows for storing up to 16384 data in a buffer, 
which can be synchronously acquired using an external 
clock. Thus we use an external DAQ card (MCC 1608FS, 
16 bit resolution, with USB interface) whose ADC opera- 
tion triggers the lock-in data storage. Such DAQ is used 
to sample reference data produced by either the pulsc- 
oximeter or an electrocardiographic connection. 



FIG. 4: a) Single input acquisition: laser frequency is tuned 
to the center of the CPT resonance, b) Differential input ac- 
quisition: CPT signals on the two arms of the sensor are sep- 
arated by more than the CPT linewidth. c) Differential input 
acquisition: CPT signals overlapping is optimized, d) Out of 
Doppler absorption resonance noise in differential input mode, 
e) Out of Doppler absorption resonance noise in single input 
mode. In the inset is reported the total rms noise N inte- 
grated betweeen 1 Hz and 30 Hz (cardio-signal bandwidth) 
relative to each trace. 

In the case of the pulse-oximeter, we arranged an IR 
led and a phototransistor in such a way to sense the 
transparency variation of a finger, and give several hun- 
dred mV signal peak-to-peak. Alternatively, the ECG 
was performed by acquiring digitally the electric signal 
collected by two electrodes placed in the vicinity of the 
heart, after a passive low-pass filtering, necessary to re- 
duce the noise aliasing. Specifically, as we use typically a 
128 S/s acquisition rate, a 18 dB/oct filter made with a 
three-stage RC cutting at 40 Hz, produced a clean signal 
about 1 mV at the QRS complex. 

In both cases, the reference signal needs some prelim- 
inary (off-line) numerical conditioning before being used 
as a trigger. Namely, some linear filtering is performed to 
remove slow drifts and high frequency noise (a 3rd order 
bandpass numerical filter from 0.2 Hz to 30 Hz is usu- 
ally suitable) . After this, the envelope of the reference is 
evaluated, and used to normalize the peak-to-peak am- 
plitude all along the registration. Finally, the normal- 
ized reference is used to produce a trigger signal which is 
adjustable in slope, threshold and hold-off, similarly to 
standard oscilloscopes. Additionally, and in particular 
when using pulse-oximeter reference, it is important to 
adjust a (negative) trigger delay, of some tenths of sec- 
ond, due to the relevant lag between the heart pulse and 
the change on the finger transparency. Both references 
showed to be suitable for data average triggering. 

The magnetocardiographic data downloaded from the 
lock-in data buffer are firstly numerically filtered with 
the aim of removing specific spectral components (such 
as 50 Hz noise) and other spectral peaks due to artifacts 
(such as mechanical vibration of the Helmholtz coils) 
then split in traces corresponding to single heart pulses, 
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from which the average pulse can be reconstructed by 
averaging. 

The averaging process is not straightforward, due to 
the fact that the pulse duration is not stable. We tried 
two different approaches to superpose and average the 
single-pulse traces. In the first one, we strengthen the 
arrays containing the single pulse traces to a given size by 
linear interpolation, to make all of them equal in length; 
while in the second procedure we add zeros at the end 
of the shorter arrays to make all of them as long as the 
longest one. It turns out that the second procedure is 
more effective and correct, due to the fact that (as can be 
seen also by the ECG traces) the traces corresponding to 
single pulse cardiac activity are not similar to each other. 

Specifically, each trace starts with a P-QRS-T complex 
having a pretty stable duration, followed by a quiescent 
interval whose duration is rather variable in time, which 
is the main responsible for the non-periodicity of the sig- 
nal. As a consequence, using ECG reference, one finds 
that the averaging is the best (and the best is the S /N ob- 
tained) when triggering on the QRS complex (the highest 
peak) with a delay suitable to include the P wave of the 
same beat, while the location where the quiescent phase 
is being cut is not relevant, provided that the T wave of 
the previous beat is not included. In the case of pulse- 
oximeter as a reference the correct trigger delay must be 
found empirically due to its intrinsic additional delay. 

In Fig. [5] is reported a reconstructed magnetic heart- 
beat, obtained by placing the sensor close to the chest of 
one of the authors. The peak is obtained from a data set 
consisting in 16384 points acquired at a sampling rate of 
128 Hz. Average is thus performed on a set of about 150 
heart-beats, corresponding to 128 sec of measuring time. 
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FIG. 5: Reconstructed cardiac pulse. The signal is averaged 
over a set of about 150 cardiac pulses acquired at 128 Hz 
sampling rate. The lock-in equivalent noise band-width is 
41.7 Hz. 



VI. CONCLUSIONS 

We have demonstrated the potential of a fully opti- 
cal, CPT based, differential magnetometers for magneto- 
cardiographic applications in unshielded environment, 
using a sample of Cs vapor at room temperature in a 
compact sensor head. We discussed simple approaches 
for magnetic-field and magnetic-field-gradient compen- 
sation, as well as off-line data analysis of quasi-periodic 
signals. The noise characterization showed that a severe 
limitation to the sensitivity would be set by the back- 
ground magnetic field fluctuation, which, thanks to its 
essential uniformity can be effectively canceled in differ- 
ential measurements. Specific features of the presented 
magnetometer, compared to others used at present lie in 
the facts that it works at room temperature and can be 
fast adjusted to detect different components of the weak 
biologic magnetic field vector, as the sensor, being purely 
optical, does not include coils, which would constitute ge- 
ometrical constraints. The simplicity of the method, its 
low cost and low maintenance cost can be a crucial factor 
for its dissemination in clinical applications. 



MAGNETIC-FIELD AND GRADIENT 
COMPENSATION 

In this appendix, we report the principles of the strat- 
egy for the resonance narrowing, shown in Fig. [31 ob- 
tained by means of magnetic field and magnetic field gra- 
dient conditioning. 

In order to be able to independently compensate 
each of the three components of the magnetic field vec- 
tor, while introducing minimum inhomogeneities, it is 
straightforward to use three pairs of coils in Helmholtz 
configuration (which for squared coils means to place 
them at a distance of about 0.544 L, L being the length 
of the coil side). 

The center of the Cs vapor cell, the core of the op- 
tical sensor, is then placed at the geometrical center of 
a "cube" composed by a set of three pairs of squared 
(180 x 180 cm 2 ) coils (180 x 0.544 is just enough to make 
possible a human body to access the cube). 

Basic arguments give that not only V • B, but also 
V x B is zero in the position of the sensor, as no current 
flows there. As a direct consequence, only five out of the 
nine elements of the matrix dBi/dxj are independent, 
i.e. only five parameters have to be controlled for their 
compensation. More precisely, a complete compensation 
requires to controll two of the three diagonal elements, 
and three of the six off-diagonal elements. 

The diagonal elements can be compensated by anti- 
Helmholtz coils. This means that two additional couples 
of coils with counter-propagating currents could be used 
to compensate dB x /dx and dB y /dy (which would guar- 
antee also the compensation of dB z /dz). In principle, 
such anti-Helmholtz coils should be at a larger distance, 
in order to produce the linear term dBi/dxi and not the 
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cubic one. But for simplicity, considering the small cor- 
rections to be made, is preferable to use the already ex- 
isting Helmholtz pairs (with which, by the way, one has 
also the advantage of producing a maximum gradient, 
given the current). 

This choice leads to simultaneously compensate the x 
and y components of the field and the three diagonal 
components of the gradient, by separately controlling the 
currents in the four Helmholtz coils for B x and B y com- 
pensation. Specifically, indicating with l x \ and I X 2 the 
currents in the coils around the x axis, B x is controlled 
by (ixi + 1x2) and dB x /dx is controlled by (I xl - I x2 ) 
(similarly for the y direction). If only B x and B y are 
to be compensated, so that only B z inhomogeneities are 
relevant a simpler control can be made with three cur- 
rents only, using pure Helmholtz configuration in x and 
y directions ((I x i = I X 2), {lyi — -^2)), and pure anti- 
Helmholtz (I z i = — / z a) configuration in the z direction. 
In our present set-up, we opted for this latter, simpler 
choice. 

The three independent off-diagonal elements of the gra- 
dient (dBy/dx, dB z /dx, dB z /dy) can be controlled with 
couples of magnetic dipoles symmetrically located far 
away from the sensor in order to produce a vanishing, 



quadrupole field in the region of the cell. Again, pro- 
vided that dBy/dx produces negligible effects because 
not responsible for B z inhomogeneities, only two of these 
three couples needs to be effectively adjusted. We actu- 
ally use dipoles oriented in z direction and located in the 
xy plane. As the dipole field gradient decreases with the 
fourth power of the distance, fixing the dipoles on the 
frame of the large Helmholtz coils would make necessary 
to use pretty large current and pretty heavy coils. We 
simplified our task by using several couples of Nd per- 
manent magnets in order to coarsely compensate the off- 
diagonal gradient components, and smaller electromag- 
nets with few hundreds mA current for fine adjustments. 
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